Transport through single-level quantum dot in a magnetic field 
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We study the effect of an external magnetic field on the transport properties of a quantum 
dot using a recently developed extension of the functional renormalization group approach to non- 
equilibrium situations. We discuss in particular the interplay and competition of the different energy 
scales of the dot and the magnetic field on the stationary non-equilibrium current and conductance. 
As rather interesting behavior we find a switching behavior of the magnetic field for intermediate 
correlations and bias voltage. 
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I. INTRODUCTION 

The investigation of transport through mesoscopic sys- 
tems has developed into a very active research field in 
condensed matter during the past decade due to their 
possible relevance for next-generation electronic devices 
and quantum computing.— The advance in preparation 
and nano-structuring of layered semiconductors^ or the 
handling of molecules respectively nano-tubes has lead to 
an increasing amount of experimental knowledge about 
such systems >^ 

The simplest realization of a mesoscopic system is the 
quantum dot."^'^'^ It can be viewed as artificial atom cou- 
pled to an external bath, whose properties can be pre- 
cisely manipulated over a wide ranged The transport 
properties of quantum dots in the linear response regime 
are very well understood from the experimental as well as 
theoretical point of view4 On the other hand, a reliable 
theoretical description of even the stationary transport 
in non-equilibrium is still a considerable challenge. 

In the present work we study the influence of an exter- 
nal magnetic field B on the stationary transport proper- 
ties of a single-level quantum dot subject to a bias voltage 
Vb aX T — Q. Quantum dots in external magnetic field 
have been the subject of interest for some time and ex- 
perimental studies of these systems have been performed 
by several groups i^i^ii^ From a theoretical point of view, 
non-equilibrium properties in magnetic field were investi- 
gated by Meir and Wingreenii combining different meth- 
ods such as noncrossing approximation (NCA) , equations 
of motion (EOM) and variational wavelenght approach in 
the limit in which the Coulomb repulsion U is very large. 
Rosch et alii^ used a perturbative renormalization group, 
which permits the description of the transport properties 
in single-level quantum dots for large bias and magnetic 
fields. Konig et al3 studied tunneling through a single- 
level quantum dot in the presence of strong Coulomb 
repulsion beyond the perturbative regime by means of a 
real-time diagrammatic formulation. However, a theory 
that allows to access intermediate coupling, bias voltage 
and magnetic field strcnghts on a unique footing is miss- 
ing so far. 

Here, we use as theoretical approach the func- 



tional renormalization group (fRG)^ extended to non- 
equilibrium, which we have introduced in a recent 
publication.^^ We discuss how the current J and the dif- 
ferential conductance G are affected by B and by the 
competition between magnetic field and bias voltage. We 
show in particular that B is responsible for a switching 
behaviour in J as function of Vb ■ Interesting are also the 
individual contributions of spin up and down electrons, 
split by the presence of the magnetic field, to the trans- 
port parameters. To test our non-equilibrium fRG we 
have moreover studied, as limiting case, the equilibrium 
situation Vb = 0, in order to compare the results of the 
imaginary-time fRG, which has been shown to provide 
a very good description of transport through mesoscopic 
systems in the linear response regimcpJ^ 

The paper is organized as follows: The next section 
is dedicated to a brief description of the model we use, 
namely the single impurity Anderson model (SIAM). 
Based on the derivation in Gezzi et al. ^'^ we will pro- 
vide an expression for the non-equilibrium fRG equations 
studied in this work. In Sec. HI we present our results for 
the stationary transport properties. The limit ^ is 
investigated first to make contact to previous work and 
find the regime, where our method is applicable. We then 
discuss how the transport parameters current J and con- 
ductance G behave as functions of an applied bias- voltage 
with and without magnetic field. Finally we consider the 
range of applicability of the non-equilibrium fRG in the 
presence oi B. A summary and conclusions will finish 
the paper. 



II. FLOW EQUATIONS FOR SINGLE-LEVEL 
QUANTUM DOT IN MAGNETIC FIELD 

We consider in the following a quantum dot consisting 
of a single level with spin quantum number cr = ±1, 
coupled to left (L) and right (i?) leads. The electronic 
states in the leads are described by a continuum of single- 
particle levels with dispersion Sj:^^, where k denotes the 
wave vector and a — L, R. Furthermore, we assume that 
the leads are always in equilibrium. The dot and the leads 
are coupled through an energy and spin independent hy- 
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bridization Va- Finally, if two electrons occupy the dot, 
they experience a Coulomb repulsion U. This situation 
is described by the long-known single impurity Anderson 
model, given by the Hamiltonianii 



H 
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The left and right reservoirs can have different chemical 
potentials through an applied bias voltage Vb = ij-l — 
fiR. An external magnetic field is taken into account by 
the Zeeman splitting of the dot level, i.e. Sa- = VaicrB /2, 
where we introduced the gate voltage Va, which controls 
the filling on the dot. Since the magnetic fields applied 
are much smaller than the Fermi energy of the leads, its 
effect on the electrons in the leads can be ignored for the 
present purpose. 

The detailed derivation of the fRG flow equations in a 
general non-equilibrium situation has been already dis- 
cussed in Ref. Since one obtains an infinite hi- 
erarchy of differential equations for the irreducible n- 
particle vertices depending on a cutoff parameter A in 
this approach jii^ii^ii^ a truncation is necessary for ac- 
tual calculations, which we realize by setting the three- 
particle vertex = 0. In addition to this truncation of 
the hierarchy of differential equations, we further neglect 
the energy dependence of the vertex function 7^ = 7^, 
which results in an energy independent single particle 
self-energy E'^. Note that in the non-equilibrium ap- 
proach, all quantities become tensors with respect to the 
branches -I- and — of the Keldysh contour. ^^'^^ In par- 
ticular the Green function and self-energy are matrices. 



which we denote by G and E, respectively. 

Last but not least, we have to specify how the cutoff 
A is introduced. As usuali^ we choose a 6-cutoff on the 
level of the non-interacting dot here, i.e. 

GloH = e(A-|c.|)G,,,o 

with Gdafli'-^) the Green function matrix of the dot for 
U — 0. Since in this case we have to deal with a non- 
interacting system, the corresponding expressions can be 
derived straightforwardly. For simplicity we assume that 
the dispersions Sj:^^ and hybridizations Va are identical 
for a = L and a = R. In that case the coupling between 
dot and leads is characterized by the quantities 
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with A'f the local density of states of the leads at the dot 
site. The result for Gdafii^^) then readsi^: 
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where /q(±lj) :— f (±(a; — /ic)) arc the Fermi functions 
of the leads. 

With these definitions and approximations, the re- 
sulting system of differential equations for the spin- 
dependent single particle self-energy and two-particle 
vertex we are going to integrate is given by (for details 
see e.g. App. B in Ref. M) 
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In expressions ([7]) and ^ 
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and the upper small greek indices refer to the branches 
of the Keldysh contour. The initial conditions at A = cx3 
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are are given by 

j.A=oo ^ Q 

and 

All other components of 7"^=°° = 0. The integration of 
the equations (O and ([8]) has to be done until A = is 
reached. 

Compared to the system obtained in Ref. 15 without 
external magnetic field the set of eqs. ([7]) and ^ show 



a more complicated structure which manifests itself in 
a spin-dependent flow for the selfenergy and the vertex. 
As Vb ^ 0, it can be showni^ that the vertex 72 can be 
parametrized as 



with spin-independent interaction parameters U . It is 
thus tempting to use the same parametrization out of 
equilibrium, too, resulting in the simpler set of equations 
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By investigating the structure of the vertex for the full 
spin-dependent flow at A = it turns out, however, that 
the parametrized version breaks exact symmetries for ex- 
changing indices. Due to the significantly reduced numer- 
ical effort for integrating Eqs. ([9]) and pO)) . it remains 
interesting to examine its value as an approximation. 

After integration of system ([7]) , ^ or ([5]) , (fTU)) we in- 
sert the resulting selfenergy into the Meir-Wingreen for- 
mula for the currenti^i^ 



J. 



/«wi(Girw-Gj/(«)) (11) 



where Gda denotes the full spin dependent one-particle 
impurity Green function. 

The conductance finally is obtained from (fTTj) by nu- 
merical differentiation with respect to Vb- 



III. RESULTS 

Together with the initial conditions for self-energy and 
two-particle vertex we can now integrate the differen- 
tial equations using a standard Runge-Kutta solver. Al- 
though the systems of equations are valid for T > 0, too. 



we restrict our discussion to the case T ~ for simplicity, 
noting that preliminary results for T > in the weak- 
coupling regime are in agreement with previous studies 
and, in the linear response regime, NRG results. 



A. Test of the parametrized flow 

Let us start by comparing the full flow according to 
Eqs. ([7]) and ([5]) with the parametrized ones Eqs. Q 
and pH]) . We observe a good agreement at small Vb for 
E ^''^ (see Fig. [T})), while deviations appear in E^^''^ 
(see Fig. [Iji). However, since |E"+'^| > |E •'^|, this 
does not affect the behaviour of experimentally relevant 
quantities even for larger B, as can be seen for the case 
of the conductance G as function of gate voltage in Fig. 
[2l Increasing the bias voltage, the agreement is still good 
for small B (see Fig. [31 full curve and circles). As soon 
as we increase the magnetic field, ImS too, shows 
deviations between the full and parametrized flow (see 
Fig. |4]). In particular, the flow of the parametrized sys- 
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at 



tem gives rise to a small real component of E+' 
A = (Fig. 13)3), which actually should not exist, and 
leads to unphysical breaking of the particle-hole sym- 
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FIG. 1: (color online) (a) Flow of real (full curve and stars) 
and imaginary part (dashed curve and crosses) of sp-'^/r, 
for Vg/T = 0, U/T = 5, Vb/T = 1, B/T = 0.116. The 
continuous lines represent the solutions of Eqs. ([71 and ([HI, 
the symbols the parametrized version Eqs. ([9]) and (|10|l . (b) 
Flow of Er+'^^/F. 




FIG. 2: (color online) Conductance G normalized to Go = 
jh as function of Vg for (7/r = 5, Vg/F = 1 and several 
values of B. 
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FIG. 3: (color online) The same parameters as in Fig.[2]except 
for VbIV = 3. 
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FIG. 4: (color online) The same parameters as in Fig.[T]except 
for Vb/F = 3 and B jV = 1. In (b) the real part of E-+ has 
been rescaled by a factor ten to make it visible. 



metry in physical quantities like the conductance, as is 
visible from the dashed and dot-dashed curves in Fig. [31 
We thus can conclude that the parametrized set of flow 
equations is a reasonable approximation at least for cal- 
culating the conductance as long as Vb and the magnetic 
field B become not too large. 



B. Equilibrium case 

Before studying the influence of a magnetic field on 
stationary non-equilibrium transport at T = 0, we first 
present results for the linear response regime, i.e. the 
limit Vb — [II — [i-R 0. This limit can serve as a test 
for the non-equilibrium fRG, because it should reproduce 
the imaginary-time fRG resultsj^^ 

In Fig. [51 we compare the contributions to the con- 
ductance from the individual spin channels calculated 
with the non-equilibrium fRG (symbols) as Vb ^ for 
different values of B to the results of the imaginary-time 
fRG obtained by Karrasch et al.— (lines). We note, that 
for Vb/F = our calculations perfectly reproduce the 
results of Karrasch et al.— . We see that, as soon as the 
magnetic field increases, starts to split into two peaks 
which reflects the field dependent shifts of the spectral 
functions for up and down spin. 

We would like to mention here that for too large U 
and magnetic field a crossing of solutions of the differ- 
ential equations appears, and the numerical procedure 
picks the one on the wrong Riemannian sheet. At present 
we do not know how to avoid this numerical instabilty, 
which limits the applicability of the method to interme- 
diate coupling parameters. We will come back to this 
problem later. 

For an extensive discussion of the linear response re- 
sults with an applied magnetic field we refer the reader 
to the seminal paper by Karrasch et al.— . 
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V^/U 

FIG. 5: (color online) Conductance Ga from the individual 
spin channels in the limit Vs — » normalized to Go = 2e^//i 
as function of Vg/U for fZ/F = 5 and several values of B. 
Symbols were obtained from the non-equilibrium fRG, lines 
from the imaginary-time fRG of Karrasch et al.— . 



C. Non-equilibrium 

Switching on the bias voltage Vb we observe differ- 
ent behaviours of the conductance G as a function of 
the gate voltage Vg depending on the competition be- 
tween the voltage and the magnetic field. In Fig. [B^-d 
we present G = G| + as function of Vq for fixed 
Coulomb interaction [//F = 5 and Vb/F = 0, 1, 3, 5 
for different values of B. We observe a drastic change of 
the conductance due to the interplay of Vb and B. In 
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FIG. 6: (color online) Total conductance G normalized to 
Go = 2e^/h as function of Vc/U for U /V = 5, Vs/F = 0, 1, 
3, 5 and several values of B/V. 

particular, we see that for small bias (Vb/F = 1, Fig. 
^p) and fields B /T = 0.116, G consists of just two peaks 
separated by U . If we now increase the value of the mag- 
netic field, the conductance at small Vq initially strongly 
increases (dashed and dot-dashed curves in Fig.lHb), the 



peaks disappear and a small plateau appears. Increas- 
ing B further, the plateau disappears and we get back 
the two peaks separated by a rather deep valley and two 
shoulders whose spacing is A « [7. 

This behaviour can be explained as follows: The spec- 
tral density for each spin channel is split by Vb into two 
peaks moving the spectral weight to higher frequencies 
and decreasing it in the region Vg ~ 0,^^ Switching on 
B will lead to a splitting ~ i? of each peak, due to field 
induced shifts of the spectral functions from the individ- 
ual spin contributions. Thus, increasing B, the two outer 
peaks move away from each other while the inner ones 
get closer and even merge, enhancing the conductance 
at Vg 0. As direct consequence we observe a small 
plateau for Vg ~ (see Fig. [6j)). Further increasing 
B, these contributions will again drift apart, leading to 
a collapse of the conductance and disappearance of the 
plateau. 

Completely different is the behaviour for Ve/r = 3 
and Vs/r — 5 (Fig. [SJ: and d). We find a monotonic 
decrease of G with the magnetic field. In addition the 
field dependence is initially weaker than in Fig. [6)d. We 
interpret this behaviour in the following way: For large 
Vb, the Fermi window in e.g. Eq. (|lip will lead to an 
averaging over a large energy region. Thus structures due 
to the magnetic field at too small energies will be washed 
out, i.e. the subtle interplay between the rearrangement 
of spectral weight due to Vb and the shifts induced by 
the magnetic field B cannot be resolved any more. 



D. Current and Conductance as function of the 
applied bias 

The non-monotonic behaviour of G(Vg = 0) for mod- 
erate Vb as function of the external magnetic field B is an 
interesting feature we want to explore in somewhat more 
detail in the following. To this end we calculated the 
current J and the differential conductance G = dJ/dVB, 
at Vc/r = and U/T — 5 a,s function of the applied 
bias for different values of B. The results are collected 
in Fig. [T] Compared to the case B = (full line) we 
see that a finite magnetic field B basically induces two 
features. First, as is also the case for B = at larger U, 
we observe a shoulder in the current, resulting in a small 
region of negative differential conductance. More inter- 
esting, however, is the appearance of an almost unitary 
conductance peak at Vb ~ B for intermediate magnetic 
fields. For large fields the features is suppressed again. 
This behaviour can be explained by noting that when 
Vb ~ B, the electrochemical potentials fJ,L,R of the leads 
are close to the split dot levels, respectively, therefore the 
tunnel probability from the leads to the dot is enhanced. 
For Vb ^ B the curves converge again to the same val- 
ues, i.e. B does not influence the current any longer, the 
behavior is dominated by the now rather wide Fermi win- 
dow. 

Thus we observe that B can act as a switch, suppress- 
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FIG. 7: (color online) Total current normalized to Jo = Go^ 
and conductance as function of the bias voltage Vb /F for 
Vc/r = 0, U/r = 5 and several values of B. 



Vs/r = 2.. 3 the conductance now overshoots the unitary 
Umit in the range 1 < B /T < 2. 

While these results are at first sight very disturbing, 
they can be traced back to a crossing of solutions of the 
differential equations in this parameter regime. As al- 
ready mentioned earlier and discussed in detail in Ref. 
[TBI , the Runge-Kutta algorithm picks the wrong solution, 
leading to the observed pathologies in the physical prop- 
erties. At very large bias Vb > U the solution of the 
differential equations does not show this problem any 
longer. 

This means that our approach is due to numerical rea- 
sons not reliable in this parameter range and possibly for 
larger magnetic fields. We are presently not aware of any 
algorithm that can avoid this numerical difficulty related 
to the appearance of a pole in the physical Riemannian 
sheet. 



ing the current for small Vb and leading to a steep in- 
crease of the current in the voltage range Vb ^ B. 



E. Range of applicability of the fRG 

Let us finally discuss in which range of parameters the 
non-equilibrium fRG furnishes reliable results. In Fig. 
[5] we show the transport parameters plotted as function 
of the magnetic field for different bias voltages in the 
weak < 1, Fig. [5^) and intermediate coupling regime 
> 1, Fig.ib). While for 1^ < 1 the fRG always con- 
verged smoothly, we see that for ^ > 1 the curves show 
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FIG. 8: (color online) Total conductance normalized to Go ~ 
2e^//i plotted versus B /F for Vg/F = and several values of 
Vb- (a) weak coupling regime C//F = 1. (b) strong coupling 
regime U/T = 5. 

a discontinuity in G (see Fig. [5^b for Vg/F = 1). In- 
creasing Vb further this discontinuity disappears but for 



IV. SUMMARY AND OUTLOOK 

In this work we applied the non-equilibrium fRG ap- 
proach to the stationary transport through a single-level 
quantum dot subject to a constant bias Vb and magnetic 
field B at T = 0. Besides the truncation of the infinite 
hierarchy of the fRG differential equations at the two- 
particle level we also neglected the energy dependence of 
the two-particle vertex function 72, resulting in an energy 
independent single-particle self energy. Guided by the 
structure of the fRG equations in equilibrium, we stud- 
ied as further possible approximation the neglect of the 
spin-dependence in the vertex function 72. This latter 
approximation further significantly reduces the complex- 
ity of the system of differential equations. We compared 
the results of this approximation with those from the 
equations maintaining the full spin dependence, finding 
that for not too large Vb and field B it actually can be 
used to accurately calculate transport properties. 

Although the truncation and especially the neglect of 
the energy dependence must be viewed as rather crude 
approximations, we have shown that one can obtain rea- 
sonable results for the transport parameters J and G 
up to the intermediate coupling regime for an extended 
regime of the parameters Vb and B. A particularly in- 
teresting observation is the switching behaviour found in 
the current for intermediate values of Vb and magnetic 
field B, which we could explain with the interplay of the 
different structures we expect in the single-particle spec- 
tra as function of the gate voltage and B. We also showed 
that for bias voltage Vb ^ we recover, as expected, the 
linear response results by Karrasch et al.— . 

As an important future step we view the introduction 
of the energy dependence in the flow equations of the 
non-equilibrium fRG. Besides curing certain deficiencies 
not discussed here (see Gezzi et al.— ), we expect that 
this step will extend the range of applicability of our 
formalism to larger magnetic fields and Coulomb inter- 
action. Furthermore, coping with the full energy depen- 
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dencc of the vertex will eventually enable us to study 
time-dependent phenomena off equilibrium, which will 
open a wide and interesting field of physical phenomena 
and applications. 
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